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SUMMARY 

An  approach  is  proposed  that  uses  coordinate  measurements  of  the  real  sur¬ 
face  of  spiral  bevel  gears  to  determine  the  actual  machine  tool  settings 
applied  during  the  gear  manufacturing  process.  The  deviations  of  the  real  sur¬ 
face  from  the  theoretical  one  are  also  determined.  Adjustments  are  then 
applied  by  machine  tool  corrections  to  minimize  these  surface  deviations. 

This  is  accomplished  by  representing  the  real  surface  analytically  in  the  same 
Gaussian  coordinates  as  the  theoretical  surface.  -  .  >•  , 
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1.  INTRODUCTION  "  ' 

The  development  of  computer-controlled  machines  has  opened  new  opportuni¬ 
ties  for  high-precision  generation  of  double-curved  surfaces  -  the  surfaces  of 
gear  teeth,  rotors,  propellers,  screws,  etc.  However,  these  opportunities  can 
only  be  realized  if  the  surface  generation  is  complemented  with  coordinate 
measurements  of  the  manufactured  surfaces.  Only  with  such  measurements  can 
the  deviations  of  the  real  surface  from  the  theoretical  one  be  determined  and 
then  minimized  by  correcting  the  applied  machine  tool  settings. 

The  Gleason  Works  (USA),  Oerlikon  (Switzerland),  Caterpillar  (USA),  and 
the  Ingersoll  Milling  Machine  Company  (USA)  are  pioneers  in  the  development  of 
computer-controlled  machines  for  generating  spiral  bevel,  hypoid,  spur,  and 
helical  gears.  Engineers  at  the  Gleason  Works  have  developed  a  system  for 
automatically  evaluating  real  gear  tooth  surfaces  that  is  based  on  measure¬ 
ments  taken  with  the  Zeiss  machine  (ref.  !).  Engineers  at  Caterpillar  have 
developed  their  own  machine  for  coordinate  measurement  and  have  used  it  for 
evaluating  and  correcting  real  gear  tooth  surfaces  (ref.  2).  Coordinate  meas¬ 
urement  of  complicated  surfaces  is  expected  to  find  wide  application  in 
industry. 
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This  paper  covers  the  following  topics: 

(1)  Numerical  determination  of  machine  tool  settings  for  a  real  surface. 
Here  it  is  assumed  that  the  deviations  of  the  real  surface  from  the  theoretical 
one  are  caused  only  by  machine  tool  setting  errors.  The  proposed  approach 
allows  the  required  corrections  of  machine  tool  settings  to  be  determined  from 
coordinate  measurements. 

(2)  Determination  of  corrections  of  machine  tool  settings  for  a  real  sur¬ 
face  with  irregular  deviations.  Deviations  can  be  caused  by  heat  treatment 
and  deflections  during  manufacturing.  The  proposed  approach  assumes  that  the 
manufacturing  process  provides  repeatable  surface  deviations  and  allows  the 
deviations  to  be  minimized  by  appropriate  corrections  to  the  machine  tool 
settings. 

(3)  Analytical  representation  of  the  real  surface  including  the  deviations 
that  remain  after  correction  by  (2).  The  proposed  approach  allows  the  real 
surface  to  be  represented  in  the  same  Gaussian  coordinates  as  the  theoretical 
surface  so  that  computer-aided  simulation  of  meshing  and  contact  of  the  inter¬ 
acting  surfaces  (e.g.,  gear  tooth  surfaces)  can  be  simplified. 

The  solution  to  these  problems  is  illustrated  by  a  numerical  example. 


2.  REPRESENTATION  OF  A  THEORETICAL  SURFACE 

A  theoretical  gear  tooth  surface  is  an  envelope  of  the  family  of  tool  sur¬ 
faces.  Methods  for  their  analytical  representation  are  well  known  and  have 
been  described  in  reference  3. 

The  theoretical  surface  Et  may  be  represented  in  a  parametric  form  in  a 
coordinate  system  St  rigidly  connected  to  Et  as  follows: 

2  9rt  9rt 

rt<u,e;q^)  €  C  <j  -  1,2,  .  .  .,n>;  u,0  G  E;  — x  — ~  ^  0  (1) 

The  designation  C^  means  that  the  vector  function  has  continuous  derivatives 
for  all  arguments  at  least  to  the  first  and  second  orders.  The  Gaussian  coor¬ 
dinates  are  designated  by  u  and  9,  and  E  is  the  area  of  u  and  0.  The 
inequality  in  equation  (1)  indicates  that  Et  is  a  regular  surface.  The 
machine  tool  settings  are  designated  by  constants  q j < j  =  1,2,  .  .  .  ,n). 

This  approach  requires  a  parametric  representation  of  a  surface  that  is 
equidistant  from  the  theoretical  surface.  Such  a  surface  is  represented  by 

r^(u,9)  +  Xn^.(u,9)  (X  *  0)  (2) 

Here 

N,  3r.  3r. 

nt(u’9)  =  ThJ'  Nt  =  “§U  x  "39  *  0 
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where  is  the  vector  of  the  surface  normal;  is  the  unit  surface  nor¬ 

mal;  and  X  is  the  scalar  that  determines  the  distance  between  the  two 
surfaces . 


3.  PRINCIPLES  OF  COORDINATE  MEASUREMENT 

A  coordinate  measurement  machine  is  supplied  with  a  probe  that  can  per¬ 
form  translational  motions  in  three  mutually  perpendicular  directions  during 
the  measurement  process.  The  probe  tip  represents  a  changeable  sphere  whose 
diameter  can  be  chosen  from  a  wide  range.  Henceforth,  we  will  consider  that  a 
coordinate  system  Sm(Xm,Ym,Zm)  is  rigidly  connected  to  ^he  coordinate  measure¬ 
ment  machine,  where  Zm  corresponds  to  the  axis  of  the  gear  (fig  u  The 
axis  of  the  probe  may  be  installed  parallel  to  Zm  (fig.  1(a))  or  perpendicular 
to  Zm  (fig.  1(b))  as  is  appropriate  -  depending  on  the  pitch  cone  angle  of  a 
hypoid  or  spiral  bevel  gear,  for  example.  The  back  face  of  the  gear  is  its 
base  plane,  and  the  origin  of  the  coordinate  system  Sm  is  located  in  the  base 
plane  or  is  related  to  it.  A  coordinate  system  S^(  x^  ,yt  's  rigidly  con¬ 
nected  to  the  gear  being  measured.  In  some  cases  we  may  assume  that  the  origin 
Ot  coincides  with  0m.  In  the  most  general  case  the  orientation  and  location 
of  St  with  respect  to  Sm  are  determined  with  two  parameters,  6  and  2. 

(fig.  2).  These  parameters  can  be  determined  by  using  the  computational  proce¬ 
dure  described  in  section  4. 

The  coordinate  measurement  machine  is  provided  with  a  rotary  table.  The 
table  allows  the  gear  to  be  installed  in  an  initial  position  with  respect  to 
the  probe.  The  measurement  data  provide  the  coordinates  of  the  center  of  the 
probe  tip  sphere. 

The  coordinate  measurement  machine  can  be  calibrated  for  a  chosen  probe 
tip  sphere  by  using  a  calibration  ring  (fig.  3).  The  initial  coordinates  of 
the  center  of  the  tip  sphere  are 


'x<0)y(0)z(0> 

.  m  m  m 

where  R  is  the  radius  of  the  calibration  ring,  a  is  the  radius  of  the 
sphere,  and  f  is  obtained  by  independent  measurement.  At  the  initial  posi¬ 
tion  the  probe  sphere  is  in  contact  with  the  calibration  ring.  The  Ym  =  0 
alignment  is  achieved  by  finding  the  Ym  position  where  equal  displacements 
±AY  of  the  probe  result  in  equal  Xm  direction  displacements.  Since  the 
probe  performs  measurements  by  translational  motion,  its  displacements  in  the 
Xm,  Ym>  ancl  axis  directions  represent  displacements  of  the  sphere  center 
from  the  initial  position. 


a,  0,  f 


(3) 


4.  DETERMINATION  OF  REAL  MACHINE  TOOL  SETTINGS 
Initial  Considerations 

All  deviations  of  the  real  surface  from  the  theoretical  one  are  assumed 
to  be  only  the  result  of  errors  in  the  applied  machine  tool  settings.  Then 
the  real  machine  tool  settings  are  determined  from  the  coordinate  measurement 
data . 
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Consider  that  the  theoretical  gear  tooth  surface  and  the  unit  surface  nor¬ 
mal  are  represented  in  coordinate  system  St  by  the  following  vector 
equations: 


u.e; 

_( 0 ) 
q1  •  • 

_  ( 0 )  ] 

-  ■ ,qn  J 

(4) 

[u.e; 

„<0> 
ql  ’  • 

•  -,qn  J 

(5) 

where  .  .  .  ,q^  represent  the  nominal  machine  tool  settings.  To  rep¬ 

resent  the  real  surface  and  its  unit  normal  in  St,  substitute  the  real 
machine  tool  settings  to  be  determined  < q l , q2 . dn^  f°r  the  nominal  val¬ 

ues  in  equations  (4)  and  (5). 


Now  consider  an  imaginary  surface  from  the  real  surface  at  a  distance 
equal  to  the  radius  of  the  probe  sphere.  This  surface  is  represented  in  St 
by  (see  eq.  (2)) 


x(e) 

xt 

=  xt(u,9; 

"j* 

+  anx(u,9; 

V 

=  A(u,0; 

V 

4e) 

=  yt(u,0; 

V 

+  any(u,0; 

V 

=  B(u ,©; 

V 

(6) 

4e> 

=  zt(u,6; 

V 

+  anz(u,6; 

V 

=  C( u ,0; 

V 

where  a  is  the  radius  of  the  probe  sphere;  A,  B,  and  C  represent  the  result¬ 
ing  functions;  and  q-j(j  =  1,  .  .  .,n)  are  the  unknown  r^al  machine  tool 
settings. 


Basic  Equations 


The  real  machine  tool  settings  are  determined  as  follows: 


Step  1 .  -  The  coordinate  transformation  from  St  to  Sm  is  based  on 
the  matrix  equation 


M  -  [Mmt]  ft] 

Here  (see  fig.  2) 

cos  6  sin  6  0  O' 

-sin  S  cos  6  0  0 

0  0  1  a 

0  0  0  1. 


CMmt] 


(7) 


(8) 


Considering  that  the  measured  coordinates  of  the  probe  sphere  center  (Xm,Ym,Zm) 
coincide  with  coordinates  on  the  equidistant  surface  represented  in  Sm, 
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'x(e)y(e) ,(e) 
.  m  m 


]  -[■ 


X  Y  Z 

m  m  mj 


Equations  (6)  to  (9)  yield 

Xm  =  A(u,9;  q.)cos  8  +  B(u,0;  q . ) s i n  8 

J  J 

Y  =  -A(u,6;  q.)sin  8  +  B(u,9;  q.)cos  8  ( 1 0 > 

J  J 

Zm  =  C(u ,0;  q j )  +  a 

Step  2.  -  The  goal  is  to  derive  equations  that  are  invariant  with  respect 
to  the  parameters  8  and  4.  The  drawing  of  figure  2  yields 

P2  =  Xm  +  Ym  =  xt  +  yt  =  a2(u-6;  V  +  B2(u,e;  q  )  (111 


Using  equation  (11)  yields 


Equations  (10)  yield 


X„  +  Y^  =  A2(u ,0;  q.)  +  B2(u,9;  q.) 
m  m  j  j 


,  A(A  -  X  )  +  B(B  -  Y  ) 

.  6_  _ m _ m_ 

Idn  2  =  BX  -  AY 

m  m 


It  is  also  evident  that 


a  =  Zm  -  C(u,0;  q.) 


Step  3.  -  Henceforth,  the  subscript  m,  indicating  that  the  coordinates 
of  a  point  are  represented  in  coordinate  system  Sm,  is  dropped.  The  number 
of  measurement  points  is  designated  by  d  and  the  index  of  a  measured  point 
by  the  subscript  p.  The  following  system  of  equations,  based  on  equations 
(12)  to  (14),  is  used  for  determining  the  real  machine  tool  settings: 

Xp  +  yp  =  A2(up,0p;  qJ  +  B2(up,ep;  q.)  (p  =  1,2,  .  .  ,,d)  (15) 


X_)  +  B  (B  - 


A  ,  (A  . 
P+1  P+ 


Xn  .)  +  B„  , (B 
P+1  P+1 


B_  , X  ,  -  A  ,  Y  . 
p+1  p+1  p+1  p+1 


(1  <  p  <  d  -  1 ) 


’P+1  P 


zn  =  C(uP+i.ep+r  V  '  C(u"’e":  n  <  P  <  d  -  i)  (17) 


p  p  j 
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From  the  results  of  measurements  for  d  points  on  the  surface,  (3d  -  2) 
equations  can  be  obtained  for  equations  (15)  to  (17)  in  2d  unknown  surface 
coordinates  (Up.Qp)  and  n  unknown  machine  tool  settings  q j ( j  =1,  .  .  ,,n). 
Thus,  determining  n  unknown  machine  tool  settings  requires 

d  =  n  +  2;  k=3d-2=3n+4  (18) 

where  d  is  the  number  of  surface  measurements  and  k  is  the  number  of  non¬ 
linear  equations  that  have  to  be  solved.  The  parameters  6  and  2,  orienting 
and  locating  coordinate  system  with  respect  to  Sm  (fig.  2),  can  be  deter¬ 
mined  from  equations  (13)  and  (14). 

If  2  =  0  (the  origin  0^  coincides  with  0m) ,  the  following  equation  may 
be  used  in  place  of  equation  (17): 

Zp  =  C(Up,9p;  qj)  (19) 

For  this  case  the  coordinate  measurements  of  d  points  on  the  real  surface 
result  in  (3d  -  1)  equations  (15),  (16),  and  (19)  in  2d  unknown  surface  coordi¬ 
nates  (up,  Op)  and  n  unknown  machine  tool  settings  qj(j  =1,2,  .  .  .,n). 

To  determine  the  n  unknown  machine  tool  settings,  use 

d  =  n  +  1  ;  k  =  3d  -  1  =  3n  +  2  (20) 


Computational  Procedure 

The  numerical  solution  of  a  large  system  of  nonlinear  equations  is  a 
complicated  problem.  For  the  case  where  2,  *  0  and  n  =  4  the  number  of  equa¬ 
tions  to  be  solved  according  to  equation  (18)  is  k  =  16.  The  system  of  non¬ 
linear  equations  can  be  solved  by  usin'g  computer  software  such  as  the  IMSL 
subroutine  DNEQNF.  However,  the  successful  application  of  this  program 
requires  a  good  first  guess  -  an  initial  set  of  unknowns  that  is  used  for  the 
first  iteration.  We  propose  a  solution  procedure  that  begins  with  a  system  of 
four  equations  using  only  the  measurements  for  two  points  on  the  surface.  This 
number  of  equations,  k  =  4,  and  the  number  of  measurements,  d  =  2,  can  be 
obtained  from  equation  (18)  considering  that  n  =  0.  This  means  that,  for  the 
first  step,  errors  in  the  machine  tool  settings  are  neglected  and  the  machine 
tool  variables  q i ,  q2,  •  •  .,qn  in  equations  (15)  to  (17)  are  set  to  the 

.  ,  ,  (0)  (0)  (0) 
nominal  values  q]  ,q2  ,  .  .  .,qn  . 


Step  1 .  -  An  initial  guess  for  the  system  of  four  equations  is  obtained 
as  follows:  (1)  an  approximate  value  for  2  is  determined  by  measurements  and 
then  (2)  neglecting  the  errors  for  machine  tool  settings,  approximate  values 
for  the  surface  coordinates  of  two  measured  points  are  determined  by  using  the 
fol lowi ng  equations : 


C(up,0p)  -  Zp  -  1  <p  - 

1  ,2) 

(21 ) 

*2<vv  *  b2<w  -  XP +  YP 

(p  =  1 ,2) 

(22) 
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Step  2.  -  Once  the  appropriate  values  of  (u,9)  for  the  two  measured 
points  are  known,  more  precise  solutions  for  the  surface  coordinates  can  be 
obtained  by  using  the  system  of  four  equations 


a2(u] ,e] ) 

+  B2(U1 ,61 )  =  X2  +  Y2 

(23) 

a2(u2,g2) 

+  B2(u2,62)  =  X2  +  Y2 

(24) 

C(u2 ,©2> 

-  C(u-|  ,9j )  =  Z2  -  Z j 

(25) 

■  Xj)  +  B1(B] 

-  Y] )  A2(A2-  X2)  +  B2(B2  -  Y2) 

(26) 

B, X.  -  A, Y, 

B0X0  -  A0Y„ 

obtained  for  equations  (15)  to  (17)  by  considering  that  d  =  2  and  neglecting 
errors  in  the  machine  tool  settings. 

Step  3.  -  The  solution  obtained  for  the  previous  step  is  then  used  as  the 
initial  guess  for  a  larger  system  of  k  =  7  equations  (15)  to  (17),  obtained 
by  considering  that  one  machine  tool  setting  is  a  variable  and  using  d  =  3 
measurement  points. 

Step  4.  -  The  number  of  machine  tool  settings  that  are  considered  as  var¬ 
iables  is  gradually  increased  until  the  exact  values  for  the  whole  set  of 
3  =  1,2,  .  .  .,n  unknown  machine  tool  settings  are  eventually  determined  by 
using  a  system  of  k  =  3n  +  4  equations  (15)  to  (17).  Knowing  the  real  values 
of  the  machine  tool  settings  allows  the  settings  to  be  corrected  and  the  devia¬ 
tions  of  the  real  surface  from  the  theoretical  one  to  be  eliminated. 

In  some  cases  the  real  tooth  surface  is  substantially  distorted  because 
of  problems  other  than  errors  in  the  machine  tool  settings.  The  procedure 
described  in  section  5  can  be  used  to  improve  the  precision  of  the  generated 
surface . 


5.  MINIMIZATION  OF  DEVIATIONS  OF  THE  REAL  SURFACE 

Now  consider  the  case  where  the  deviations  of  the  real  surface  from  the 
theoretical  one  are  caused  by  many  factors  -  not  just  errors  in  machine  tool 
settings.  It  is  assumed  that  the  process  of  manufacturing,  including  heat 
treatment,  provides  repeatable  deviations.  The  surface  deviations  may  then  be 
compensated  for  (but  not  made  zero)  by  using  directed  corrections  to  the 
initially  applied  machine  tool  settings.  The  procedure  for  determining  the 
corrected  machine  tool  settings  is  based  on  minimizing  an  objective  function 
in  n  variables  that  describes  the  real  surface  deviations;  n  is  the  number 
of  machine  tool  settings  to  be  corrected. 

The  following  stages  of  solution  to  this  problem  are  considered:  (1) 
determination  of  the  orientation  of  the  coordinate  system  with  respect  to 
Sm  (fig.  2);  (2)  determination  of  the  deviations  of  the  real  surface;  and  (3) 
derivation  and  minimization  of  an  objective  function  (described  later  in  this 
paper) . 
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Determination  of  Orientation  Parameter  6 


The  gear  is  assumed  to  be  installed  on  the  coordinate  measurement  machine 
flush  against  a  base  plate  such  that  the  parameter  2,  is  known  (for  conven¬ 
ience,  take  2  =  0).  Furthermore,  it  is  assumed  that  the  probe  can  be  ini¬ 
tially  installed  at  a  specified  position  and  that  the  gear  can  be  rotated  until 
it  contacts  the  probe  at  this  position.  The  specified  position  of  the  probe  is 
determined  from  the  equations  for  tangency  of  the  probe  with  the  theoretical 
surface  at  a  selected  surface  reference  point  M.  This  point  can  be  chosen  as 
the  middle  point  of  the  surface  and  is  given  by  surface  coordinates  (u*,0*) 
that  determine  M  and  the  unit  surface  normal  at  M. 

The  equations  for  tangency  of  the  probe  sphere  and  the  theoretical  surface 
at  M  are  derived  by  using  equations  (6)  to  (8)  and  considering  that  the 
parameter  2  in  matrix  (8)  is  zero.  Equations  (6)  to  (8)  represent  a  system 
of  three  equations  in  four  unknowns:  Xm,  Ym,  Zm,  and  S.  If  Ym  =  0,  then 
Xm,  Zm,  and  S  can  be  determined  from  the  following  set  of  equations  derived 
from  equations  (12)  to  (14): 


X2  =  A2 
m 


u*,Q*;  q^0)]  +  B2[u*,e*;  q^ 

Zra.c[U*,9*;  ,«>] 


(0) 


(27) 


tan  = 


A2[U*,e‘;  q<0)]  -  XmA 

u*,©*;  q(10) 

X  R [n* ,0* ;  q^0) 
m  l  j 

Here  (u*,  0*)  are  the  surface  coordinates  of  the  chosen  reference  point  M; 

qj0)(j  =1,2,  .  .  . ,n)  are  the  theoretical  machine  tool  settings;  and  ( Xm , 0 , Zm) 

are  the  coordinates  of  the  center  ef  the  probe  sphere  that  is  to  be  installed 
on  the  coordinate  measurement  machine.  If,  in  the  installation  of  the  probe, 
the  actual  coordinates  (Xm,0,Zm)  differ  from  the  calculated  ones,  the  values 
of  u*,  0*.  and  S  must  be  corrected  by  using  equation  (27). 


Determination  of  Deviations 


Once  the  orientation  parameter  S  is  known,  the  theoretical  surface  can 
be  represented  in  coordinate  system  Sm.  Consider  that  Rp(Xp,Yp,Zp)  is  the 
position  vector  of  the  center  C  of  the  probe  sphere  and  that  this  point  lies 
on  the  line  of  action  of  the  theoretical  surface  normal  at  a  point  T  of  the 
theoretical  surface  (fig.  4).  The  vector  equation  is 


R 


P 


u  >0  ; 
P  P 


."pK 


(0)1 
j 


(28) 


Here  ('jp,0p)  are  the  curvilinear  coordinates  of  point  T  of  I*,  qj0) 

<j  =  1,2,  .  .  ,,n)  are  the  theoretical  machine  tool  settings;  and  Xp  is  the 
distance  between  T  and  C. 
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Equation  (28)  yields  a -system  of  three  equations  in  the  unknowns 
(Up,0p,Xp)  that  can  be  represented  as 


X  -  x 

P  P 

r  .  (o)i 

a  ,  o  ;  q  • 

Lp  p  J 

Q. 

t 

Q_ 

>-  1 

-  <0)1 
LvV  "i  J 

(■  r\  (0  )' 

n  u  ,0  ;  q . 
xpL  P  P  J 

1  n  [u  ,0  ;  q(0) 

l  ypL  p*  p  hp  J 

ZP  "  2P 

r  n  (0)i 

u  ,0  q  - 
L  P  P  3 _ 

n 

2P 

(0)1 

UrV°n:  qx 

L  M  J  -1 

(29) 


where  p  =  (1,2  .  .  .,d);  j  =  (1,2  .  .  ,,n);  and  all  coordinates  are  repre¬ 
sented  in  coordinate  system  Sm. 

The  procedure  of  computation  is  as  follows: 

Step  1 .  -  The  surface  coordinates  for  a  measurement  point  can  be  deter¬ 
mined  from  the  equations 


f 


L  ypv  P 


n  Y  - 


xpv  p 


=  0 


(30) 


u  ,  0 
P  P 


[n  (V 

-  y  )  -  n 

(l  ~  0 

zpk 

P  *V 

ypv  p  p/J 

(31  ) 


Step  2.  -  The  value  of  Ap  can  be  determined  by  using  any  of  the  three 
equations  of  system  (29) . 


Step  3.  -  The  deviation  Ap  of  the  real  surface  from  the  theoretical  sur¬ 
face  can  be  determined  by  considering  that  the  point  of  tangency  of  the  probe 
sphere  with  the  real  surface  lies  at  a  distance  a  equal  to  the  radius  of  the 
probe  sphere  (fig.  4).  Thus,  the  position  vector  Mp  of  the  point  cr'  the 
real  surface  Er  is  given  by 


where 


(0) 

VV 


+  A  n 
P  P 


w«50)] 


Ad  =  Xp 


(32) 

(33) 


Image  Surface 

The  process  of  manufacturing  initially  provides  the  theoretical  surface 
Et.  but  owing  to  unknown  factors  the  surface  E{-  becomes  distorted  into  the 
real  surface  Er  represented  numerically  by  equation  (32).  If  the  deviations 
are  repeatable,  tne  inevitable  distortion  can  be  prepared  for  and  the  devia¬ 
tions  in  the  final  manufactured  surface  minimized  by  generating,  not  the  theo- 

* 

retical  surface  E^,  but  an  image  surface  Er  represented  by 
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(34) 


N  (u.0  >  =  r 
P  P  P  P 


q 


(0)’ 


L  p*  p*  Hj 


)pnp[up 


Comparing  equations  (32)  and  (34)  shows  that  Er  is  the  image  of  Er 

reflected  through  the  theoretical  surface  E^  (see  fig.  4).  Henceforth, 

* 

Er  will  be  referred  to  as  the  image  surface. 

* 

Although  Er  is  the  surface  to  be  generated,  it  cannot  be  provided 

* 

exactly  with  the  existing  generation  process.  The  surface  Er  can  only 

* 

be  approximated  with  a  corrected  surface  Ef  represented  by  the  vector 
function 


r  (u„,e  •  q.) 

P  P  P 

Here  the  designation  qj  instead  of  qj^  means 
settings  must  be  applied  for  generating  Zf.  The 
Ej-  can  be  determined  by  using  equations  that  are 


(35) 


that  new  machine  tool 

* 

deviation  of  Er  from 
similar  to  equation  (29) : 


X*  ~  xn(un-  q-i> 

P  P  P  P  J 

n  (u*,9*;q.) 
xp  p’  p’Mj 


vu;-v  v 


Z*  -  z  (u*,e*;  q.) 

P  P  P  P 

n  ( u* ,8* ;  q  . ) 
zp  p’  p’ 


A*(u*,0*;  q.) 
P  P  P  J 


(36) 


*  * 

The  designation  (up,0p)  instead  of  (up,0p)  means  that  the  surface  coordinates 

associated  with  each  measurement  p  will  be  changed,  since  new  machine  tool 
settings  are  appl ied. 


Determination  of  Corrected  Machine  Tool  Settings 

The  goal  is  to  determine  new  machine  tool  settings  such  that  the 
★  *  * 
differences  Ap  (p  =  1,2,  .  .  .,d)  between  the  surfaces  Er  and  E^  are  mini¬ 
mized.  The  solution  to  this  problem  is  based  on  the  minimization  of  the  objec¬ 
tive  function 

d 

F  a  (A*)2  (37) 

P  P 
P=1 

where'  F  is  a  function  of  n  variables,  the  machine  tool  settings 

.,n);  and  ap  are  the  weighting  coefficients.  The  use  of 
weighting  coefficients  allows  smaller  deviations  to  be  provided  at  points  where 
higher  precision  is  required. 
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Numerical  Example 


The  results  of  measuring  a  Formate  hypoid  gear  are  represented  in 
figures  5  and  6.  The  number  of  measured  points  is  d  =  45.  Figure  5  illus¬ 
trates  the  deviations  Ap ( p  =1,2,  .  .  .,45)  of  the  real  surface  from  the 
theoretical  one.  The  locations  of  measured  points  are  represented  on  a  plane, 
and  the  deviations  Ap  are  shown  as  normal  displacements  from  the  plane.  Fig¬ 
ure  6  represents  the  same  data  on  a  plot.  Each  latitudinal  cross  section  of 

figure  5  is  represented  by  a  line  segment  in  figure  6.  Figure  7  shows  the 

* 

deviations  of  the  imaginary  surface  Er  to  which  the  corrected  theoretical 

* 

surface  is  to  be  fitted. 

★  ★  * 

The  deviations  Ap  between  the  surfaces  E^  and  Er  were  minimized  num¬ 
erically  through  use  of  equation  (33)  and  the  quasi-Newton  method  (ref.  4)  as 
implemented  in  the  I  MS  L  subroutine  UMNIF.  The  weighting  coefficient  ap  =  1 

was  used.  Figure  8  shows  the  deviations  from  the  theoretical  of  both  the 

*  * 

image  Sr  and  the  fitted  image  surface  E^.  Figure  9  is  a  three-dimensional 

★ 

representation  of  the  deviation  of  the  fitted  image  surface  Et  from  the 
theoretical  one. 


6.  ANALYTICAL  REPRESENTATION  OF  THE  REAL  SURFACE 

As  was  mentioned  previously,  two  cases  of  deviation  from  the  real  surface 
may  be  considered:  (1)  when  the  deviations  are  caused  only  by  using  the  wrong 
machine  tool  settings,  and  (2)  when  the  deviations  are  caused  by  many  unknown 
factors.  In  the  first  case  the  real  surface  can  be  represented  by  the  same 
equations  as  the  theoretical  one  just  by  substituting  the  theoretical  settings 
for  the  real  ones  (section  4).  In  the  second  case  the  deviations  can  be  mini¬ 
mized  by  correcting  the  applied  machine  tool  settings,  and  the  problem  is  to 
represent  analytically  the  new  real  surface  obtained  by  manufacturing  a  new 
gear  with  the  corrected  machine  tool  settings.  The  goal  is  to  represent  this 
surface  as  the  sum  of  two  vector  functions: 


Here  the  vector  function 


~  (0) 

u,Q;  qj 


rt [u .9;  q. 


(0)i 


3 


i  s 


+  Ar( u ,0)  ( 38) 

the  same  as  that  for  the  theoretical 


surface;  and  Ar(u,6)  is  an  analytical  vector  function  of  the  deviations  of 

the  new  surface  from  the  theoretical  one.  Even  though  the  new  gear  was  not 

(0) 

manufactured  by  using  the  nominal  settings  qj  ,  its  representation  is  based 


on 


r 1 1 u , 0 ;  q^0) 


because  the  corrected  machine  tool  settings  were  designed  to 


bring  the  final  manufactured  surface  closer  to  r^ 


u  ,9; 


If  new  measurements  and  the  procedure  of  section  5  are  used,  the  new  sur¬ 
face  can  be  represented  numerically  by  equation  (32).  Since  the  numerical 


deviations  Ap  can  be  represented  analytically  as  A(u,0),  the  new  surface 
can  be  represented  by 


r(u,6>  =  rt  u,6;  qj0) 


A( u ,9)n^ 


(0) 
3 


(39) 


However,  the  equations  for  the  surface  normal  of  the  surface  represented 
earl i er  derived  from 

N 


nr  = 


Nr 


w  _  3r 

Nr  ”  3u  x 


3r 

30 


(40) 


become  too  complicated  for  practical  use  in  tooth  contact  analysis.  For  this 
reason  it  is  simpler  to  consider  the  deviations  of  the  real  surface  as  meas¬ 
ured  along  a  vector  of  constant  direction  -  for  instance,  along  the  unit  sur¬ 
face  normal  n°  to  the  theoretical  surface  at  the  reference  point.  Such 
deviations  can  be  determined  by  using  the  following  equations,  similar  to  equa¬ 
tions  (29)  and  (33): 


-  ».(vv  ^0)) 


n 


n! 


=  X  (u  .9  ) 
P  P  P 


(41 ) 


A;  -  x;  -  a  (42) 

Here  n°(nx,  ny,  nz)  is  the  unit  normal  to  the  theoretical  surface  at  the  ref¬ 
erence  point;  tXp.Yp.Zp)  are  the  current  coordinates  of  the  center  of  the 

probe  sphere;  a  is  the  radius  of  the  probe  sphere;  and  Ap  is  the  surface 
deviation  in  the  direction  of  n° . 


Given  that  the  numerical  deviations  Ap  (Up,0p)  p  =  1 ,2 . d  have 

been  determined,  the  problem  of  fitting  a  function  A°(u,9)  to  the  numerical 
data  can  be  approached  in  a  number  of  ways.  A  common  approach  to  such  a  prob¬ 
lem  is  to  fit  piecewise  polynomial  functions  to  the  data,  such  that  the  result¬ 
ing  function  goes  through  all  of  the  data  and  is  smooth  at  the  boundary  between 
pieces  (ref.  5).  This  approach  is  not  applicable  to  the  given  problem  for  two 
reasons.  First,  the  large  number  of  segments  of  polynomial  functions  further 
complicates  the  analytic  representation,  making  tooth  contact  analysis  and  sim¬ 
ulation  of  meshing  more  difficult.  Second,  it  is  unreasonable  to  assume  that 
the  data  are  so  precise  that  the  desired  analytical  representation  must  match 
the  data  at  every  point. 

A  good  solution  to  this  problem  is  one  that  satisfies  the  two  conflicting 
goals.  The  solution  must  result  in  a  simple  expression  for  A°(u,9)  and  must 
represent  the  numerical  data  with  good  accuracy.  We  thus  propose  to  determine 
A°(u,0)  by  using  linear  multiple  regression  analysis  (ref.  6).  In  particular, 
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if  a  solution  of  the  following  form  is  assumed: 

A°(u,0>  =  bQ  +  b^'u.0)  +  b2f 2( u ,9)  +  .  .  .  +  bkfk(u,0)  (43) 

where  bg,  b i ,  .  .  .,b|<  are  undetermined  coefficients  and  f i ( u ,0) , 
f2<u,9),  .  .  .  ,f|^(u,0)  are  any  set  of  linearly  independent  functions  of  u 
and  0  not  involving  unknown  parameters.  Then  linear  multiple  regression  can 
be  used  to  determine  the  coefficients  bQ,  b ] , .  .  . ,b^  that  provide  a  least- 
squares  fit  of  equation  (43)  to  the  numerical  data.  Computer  programs  for  lin¬ 
ear  multiple  regression  analysis  exist  in  various  software  packages  such  as 
the  IMC  statistics  library  and  the  SPSS  statistics  program.  This  software 
can  soi^e  equation  (43)  for  the  unknown  coefficients,  provide  statistics  on 
the  expected  error  of  the  approximation,  and  automatically  test  various  combi¬ 
nations  of  user-supplied  functions  to  allow  the  user  to  select  the  best  subset 
of  functions.  Although  the  technique  allows  for  general  functions  f](u,0), 
f2(u,0),  .  .  .  ,f|^(u,0),  real  data  suggest  that  the  deviations  can  be  suffici¬ 
ently  represented  by  a  second-order  polynomial: 

A°(u,G)  =  bQ  +  b1 u  +  b90  +  b3u2*  b4©2  +  bBu0  (44) 


7.  CONCLUSIONS 

In  this  paper  an  overview  has  been  presented  describing  the  interrelation¬ 
ships  between  gear  geometry,  manufacture,  and  measurement.  A  methodology  also 
has  been  presented  to  improve  convergence  between  theoretical  and  manufactured 
surfaces  by  adjusting  the  machine  tool  settings  during  manufacture.  This  proc¬ 
ess  can  be  carried  out  at  several  stages  of  gear  manufacture  (e.g.,  cutting 
and  grinding)  if  so  desired.  The  methodology  can  even  be  used  to  decrease  the 
distortion  effect  of  other  manufacturing  processes  such  as  heat  treatment. 

The  following  specific  results  were  obtained: 

1.  A  process  was  devised  for  determining  the  real  machine  tool  settings 
based  on  coordinate  measurement  of  the  manufactured  gear. 

2.  A  procedure  for  minimizing  deviations  of  the  real  surface  by  correct¬ 
ing  the  machine  settings  was  developed. 

3.  An  approach  for  analytically  representing  the  real  gear  surface  was 
developed . 
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